arXiv: 1505.03149vl [hep-ph] 12 May 2015 


RUNHETC-2015-03 


Thermal Goldstino Production with Low Reheating Temperatures 

Angelo Monteux 1 ’^ and Chang Sub Shin 1 ! [j] 

l New High Energy Theory Center, Department of Physics and Astronomy, 

Rutgers University, Piscataway, NJ 08854, USA 

Abstract 

We discuss thermal production of (pseudo) goldstinos, the Goldstone fermions emerging from 
(multiple) SUSY breaking sectors, when the reheating temperature is well below the superpartner 
masses. In such a case, the production during matter-dominated era induced by inflaton decay 
stage is more important than after reheating. Depending on the SUSY breaking scale, goldstinos 
are produced by freeze-in or freeze-out mechanism via 1 —> 2 decays and inverse decays. We 
solve the Boltzmann equation for the momentum distribution function of the goldstino. In the 
freeze-out case, goldstinos maintain chemical equilibrium far after they are kinetically decoupled 
from the thermal bath, and consequently goldstinos with different momentum decouple at different 
temperatures. As a result their momentum distribution function shows a peculiar shape and 
the final yield is smaller than if kinetic equilibrium was assumed. We revisit the cosmological 
implications in both A-parity conserving and i?-parity violating supersymmetric scenarios. For the 
former, thermally produced goldstinos can still be abundant enough to be dark matter at present 
times even if the reheating temperature is low, of order 1 GeV. For the latter, if the reheating 
temperature is low, of order 0.1 — 1 GeV, they are safe from the BBN constraints. 
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I. INTRODUCTION 


When considering consequences of early universe cosmology, it is customary to assume 
that all interesting phenomena, such as dark matter production, take place after the infla¬ 
tionary period, and in particular after reheating, which marks the moment when the energy 
density of the inflaton decay products (radiation) dominates over the inflaton energy den¬ 
sity. This is justifiable, as any primordial abundance was inflated away and large amounts 
of entropy are injected during the inflaton decay stage, further diluting any other particle 
produced during reheating. 

However if the reheating temperature T Rl the maximum temperature of the thermal bath 
in the radiation-dominated era, is well below the mass scale relevant for the particle pro¬ 
duction, the abundance produced during the inflaton decay stage could be more important 
than the one produced after reheating. There are several early works to calculate the pro¬ 
duction of WIMP-like particles and its phenomenological consequences during this period 
when the reheating temperature is low enough M- In this case, the inflaton decay stage 
can be approximated as matter-dominated era with constantly injected radiation, and the 
thermal production can be calculated independently from the specific inflation model. The 
higher temperature can be achieved during this era and naive (Boltzmann) exponential sup¬ 
pression fpr the abundance is replaced by power suppression. On the other hand, in the 
studies of thermal production of super-weakly interacting particles (SWIMP), e.g. graviti- 
nos and axions, such consideration have been ignored so far, with many works focusing on 
the production at high TrQ 

One well motivated class of SWIMP is a goldstone- or goldstino-like particle, (, whose 
interaction to the visible sector is suppressed by the symmetry breaking scale, and this 
scale is much higher than its mass. This class of particles is particularly interesting in 
the sense that their elastic scattering rate is doubly suppressed by the symmetry breaking 
scale compared to their production rate. Therefore kinetic decoupling happens earlier than 
chemical decoupling. Actually this does not give any difference if the production happens at 
high Tfi or during radiation dominated era, because chemical interactions also give a thermal 
distribution for ( as f((p) oc e ~ p ^. However if Tr is well below the relevant particle’s mass 
1 Ref. 4] studied the production of gravitinos taking into account a proper treatment of the reheating 
process, with Tr m. In [5], the axino freeze-in thermal production from the neutralino decays is 
shortly discussed for low reheating temperature after thermal inflation. 
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that produces £, and the production during early matter-dominated era is important, the 
situation is different. As we will show below, the production rate at low temperature is 
quite momentum dependent, and f((p) is no longer proportional to the thermal distribution. 
Their inverse decay rate becomes more complicated, so a detailed treatment of Boltzmann 
equations is needed especially when the symmetry breaking scale is small such as in low-scale 
gauge mediation. 

In this article, we study thermal production of SWIMPs at low reheating temperature 
considering all the aspects above. We will focus on the production of the goldstino in 
supersymmetric theories, but our treatment would apply equally to similar particles. In 
theories where local supersymmetry (SUSY) is broken spontaneously, the resulting goldstino 
is incorporated in the spin-1/2 degrees of freedom of the gravitino, which acquires a mass 
777 . 3/2 = F/y/?>Mp 1 F being the scale of SUSY breaking and M P = (StiGn)^ 1 ^ 2 = 2.4 x 
10 18 GeV the reduced Planck mass. Because the coupling between the goldstino and the 
visible sector is suppressed by 1 /F, their mass and thermal production rate are tightly 
related. The relation between the goldstino mass and the interaction strength can be changed 
if there are multiple sectors with independent SUSY breaking interactions: multiple goldstini 
[6] arise and the scenario deserves more phenomenological interest. In particular, while 
one goldstino is still eaten by the gravitino, there are uneaten goldstini, whose mass is not 
unambiguously set: at first, in Ref. [5], it was shown that they generally acquire a mass 2t77 3 / 2 
from supergravity interactions. This was extended in Ref. [7], where it was computed that 
the goldstino mass could vary around 2777 3 / 2 depending on the SUSY breaking dynamics, 
and in Ref. [8j, where it was shown that, even in the global SUSY limit, the uneaten 
goldstinos could receive a large mass (up to (9(100) GeV) if multiple SUSY breaking sectors 
communicate with the Standard Model (SM) via gauge interactions: in this last case, the 
SM fields take the role of messenger fields in mediating SUSY breaking between the two 
sectors, and large masses for the goldstino can be achieved. The main interesting change 
to the standard prediction is that the mass of the uneaten goldstinos can be parametrically 
larger than the gravitino mass, and can be taken as a free parameter depending on the 
specifics of the SUSY breaking dynamics. 

This article is structured as follows: in Section [TTJ we recall and discuss in more details 


the goldstino interactions. In Section III, we discuss the early matter-dominated era and 
study the production of goldstinos when the reheating temperature is well below the super- 
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partner scale, Tr < m/10. If the F-terrn is small, the goldstino is in thermal equilibrium 
with the MSSM sector during the matter-dominated era and the freeze-out temperature is 
momentum-dependent. This happens because the 2 4 2 scattering that would bring the 
goldstinos in kinetic equilibrium has already frozen out, and only the high-momentum gold- 
stinos can efficiently inverse-scatter into the MSSM thermal bath. We numerically solve 
the Boltzmann equation for the momentum distribution function (instead of the equation 
for the number density), and fold an analytical solution that reproduces well the numerical 
results. We compute the resulting goldstino number yield Yq = n^/s, which is reduced with 
respect to the results that one folds assuming kinetic equilibrium. Otherwise, for large F- 
term, the goldstino is slowly produced via superpartners decays (freeze-in). We also consider 


non-thermal production from direct inflaton decays. We continue in Section IV, where we 
consider the late-time implications of the produced goldstinos. If /2-parity is conserved, 
they can be cold dark matter for reheating temperatures of order 1 GeV, and overdose the 
universe for larger Tr; on the other hand, if /2-parity is violated, they can decay and will 
typically interfere with Big Bang Nucleosynthesis (BBN). We derive bounds on the reheating 
temperature in the range of 0.5 — 10 GeV for given RPV couplings and goldstino masses. 
We present final remarks and conclude in Section |V) 


II. GOLDSTINI INTERACTIONS 

If there is a single SUSY breaking sector, then the corresponding massless fermionic 
degree of freedom, the goldstino, forms the spin-1/2 degrees of freedom of the gravitino, 
which has a mass 1713/2 = F/y/3Mp. The situation becomes more complicated if there 
are multiple sectors with independent SUSY breaking dynamics [6]: then, in the limit in 
which each sector is decoupled from the others, each enjoys its own SUSY algebra and, 
if SUSY breaking occurs independently in each sector, multiple massless goldstinos would 
arise. Introducing gravitational interactions, the multiple SUSY algebras are broken down 
to a diagonal subgroup, and only a linear combination of those goldstini is eaten by the 
gravitino. As discussed above, the mass of the uneaten goldstinos is not proportional to the 
value of the corresponding F-ter ms, and can be taken as a free parameter. For example, 
consider the case of two SUSY breaking sectors with hierarchical F-terms Fi F 2 . The 
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gravitino mass is 


ygr+if ^ F\ 

1713/2 y/3Mp ~ v^Mp’ 


( 1 ) 


while the different goldstini interactions are suppressed by the different F-terms. If each 
sector contributes SUSY breaking scalar mass squares 2 , and gaugino masses m,u ,2 we 
have the following coupling to a matter multiplet and gauge multiplet (A a , A“): 


C 


int 


m 4>l ii: f 

: 7 ~# <t>\ ~ 


im x 1 


e// 


\/2F f 


e// 


+ h.c.. ( 2 ) 


Here we have denoted by ij the longitudinal component of the gravitino (the eaten goldstino) 
and by ( the uneaten goldstino. F e ff = \J Ff + F 2 2 ~ Fi, and F^ ~ F 2 are the corresponding 
F-terms of ?/ and £, respectively. When \rhi\ < |m 2 |, the interactions of £ are enhanced with 
respect to those of 77 , while the mass of ( can be kept as a free parameter, much greater 
than F 2 /Mp. 

In the following, we will focus on the uneaten goldstino ( and we will simply refer to it 
as the goldstino. We will denote its F-term by F<y It is understood that the spectrum also 
includes the gravitino, which can be neglected because of its suppressed interactions. 

When the temperature falls well below the sparticle mass scale, m, the effective interac¬ 
tions between the goldstino and other light particles are useful to estimate (non-resonant) 
elastic scattering rates. After integrating out sparticles, two-goldstino interactions are given 
as P fl4] 


£«// = Acync* - (3) 

If F-parity violating interactions are introduced in the superpotential as W = A, ,/,<!), 
there are single-goldstino interactions even after integrating out sfermions [15], 

C ef7 = 12 ~ 2 f + h.c.. (4) 

This interaction is not only important to determine the life-time of the goldstino, but also 
to calculate thermal production at low temperatures. 
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III. GOLDSTINO PRODUCTION IN MATTER-DOMINATED ERA 


A. Cosmology during Reheating 


The reheating temperature is usually referred to as the temperature of the thermal bath 
when the inflaton decays.^] and it is usually defined as H{Tr) = Tj, assuming that the 
radiation energy density, p r , dominates the universe. In fact, when t = 1 /Tj, where Tj is 
the decay rate of the inflaton, the inflaton energy density is still not-negligible and such a 
temperature was achieved only during matter domination. Thus entropy injection is still 
occurring after Tr. One more precise way to define Tr is to consider the asymptotic behaviors 
of p r in both matter- and radiation-dominated epochs, as shown in Fig. [lj 

30)(a/a fl )- 3 / 2 for T » Tr, 
(7r 2 9 .(T)T‘/30)(a/a a )- 4 for T < T ri , 


Pr 


30 


(5) 


where g*(T ) is the effective number of massless degrees of freedom at the temperature T, a 
is the scale factor and a# is the scale at which two asymptotic lines meet. Numerically we 
get 

90 \ 1/4 


T r = 0.7 




( 6 ) 


x 2 gJT R ) 

This will be our definition of the reheating temperature. 

On the other hand, the thermal bath in the matter-dominated era reached a higher 
temperature, 


„ ( 24 9 ,(r H ) y /! 7M 

MAX V5V 9 2 (T„.«)J frj 


1/2 

^ I Tr, 


(7) 


where p = Vp 1 is the energy scale of inflation, and slow-roll was assumed for this simplified 
expression. Because of the continuous entropy injection during the matter-dominated period, 
the temperature scales as T/Tmax oc a~ 3 ^ 8 instead of the typical scaling of the radiation- 
dominated era, T oc gU 1 (see Fig. [I]). We refer to Ref. [2] for the details of the computations 
leading to Tmax■ If is sufficient here to note that Tmax Tr can be much larger than 
2 While we will be referring to the inflaton as the field dominating the energy density of the universe in 
the early matter-dominated era, there are other cases that reproduce the same scaling of the energy 
density, px ~ a~ 3 ; a typical examples would be the late decay of moduli, which typically result in very 
low reheating temperatures. In any case, we will refer to the reheating temperature as the maximum 
temperature achieved in the last radiation-dominated era, the one leading to Big Bang Nucleosynthesis. 
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FIG. 1. Left: schematic behavior of the energy densities of the inflaton and of thermal radiation 
(MSSM fields) during the matter-dominated era; as p r ~ T 4 , we have Tmax Tr. Right: 
Evolution of the goldstino yield with fixed reheating temperature Tr = 20 GeV (corresponding 
to a R = 5430), squark mass rriq = 1 TeV and two choices of F \= (100 TeV) 2 ,(5000 TeV) 2 , 
corresponding to freeze-out and freeze-in, respectively. 


the TeV (sparticle mass) scale, for typical inflation scales /i ~ 10 8 GeV — 10 16 GeV. The 
hierarchy of relevant scales can be summarized as 

mm ( , T r <7/1 = (T(TeV) < T MAX , (8) 

where 0 are SM particles, m^ is the mass of £, which is taken as the NLSP by assuming 
7773/2 < •C m. We consider the thermal production of goldstinos at T < Tmax, when Tr 
is taken low, Tr = (P(GeV). 

The goldstino can be produced (and annihilated) via three different channels: i) scatter¬ 
ing: 0 + ++ ( + -0, A + A -H- C + A, ..., ii) (inverse) decay: 0 C + 0, A <(->■ £ + A ..., 

and iii) RPV scattering: 0 + 0 ( + 0 (here 0, 0, A, A ^ respectively stand for SM fermions, 

sfcrmions, gauginos, and gauge bosons). 

The Boltzmann equation for the number density n^ can be written as 

h c + 3 Hn ( =g ( J J^ c \k] ( 9 ) 

= ( CT ++-+b ,i; ) cftA^ n <t> nA ii. — ( a i^4 >a ij. v )^ + ''' (0 

+ (iy*^ £ij> n c n ip T • • • (ii) 

((j^ , (iii) 
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( 10 ) 


UiiPi)--- UAvj) 0 - 


where 7 = l/y/l - v 2 , f d :i p ; /(2-) : '‘ p. t ) and 

, n \ = 1 /‘ 9^Vi _ fl^ 3 P,y 

■ ■ ■ n^ j J (27t) 3 (27t) 3 

Here, is defined in the rest frame. We take n ^ = n \^ for T = (j), ip, A, because 

they are all in thermal equilibrium for the relevant temperature scale. One might rewrite the 

inverse scattering and decay terms in the the RHS of the Boltzmann equation, and Eq. (§ 

becomes 

h c + 3 Hn c = ( K {a 4>A ^v) T n*n e 2 ii + (r^ c ^ 1 ) T n eq + ...)(!- njn \: q ) 


= F prod K 


eq 


n 


c). 


(ll) 


where (• • • )r denotes thermal average. In the treatment of the Boltzmann equation, we 
have neglected quantum-statistical effects (Pauli-blocking/Bose-enhancement), as we have 
fi< 1. For low F-terrn, these (goldstino-number-changing) interactions can be in chemical 
equilibrium at high temperatures, until the interaction rate drops below the Hubble rate 


and the process freezes out. The production rate Vj trod of Eq. (11) determines the chemical 
freeze-out (decoupling) temperature, T/. 0 ., defined by 3 H(Tf Q ) = Tp rod . Depending on the 
value of T/.o., two distinct situations for the goldstino production are possible, displayed in 
Fig. [T] 


• freeze-in: for Ty. 0 . m, goldstino interactions were not in thermal equilibrium when 
superpartners were abundant; goldstino abundance is gradually increased to a maxi¬ 
mum, after which they are diluted. 


• freeze-out: for Tf n -C m, goldstinos maintain chemical equilibrium with the super¬ 
partners until the latter are not abundant. 

Since the production at high temperatures T > m is diluted away, we can only focus 
on the production for T < m. At such low temperature, the Boltzmann equation can be 
much simplified by ignoring the scattering contribution to the production of ( [HU [17] ; for 
Tf.o. m, the freeze-in production by </> —>■ £ + -0 (A —>■ £ + A M ) dominates over the diluted 
freeze-out contribution. One can also neglect the inverse decay term given by n^/n eq in the 
last expression of the first line of Eq. 0 - 

For the freeze-out case, the situation is more subtle. It should be noted that the factoriza¬ 
tion by (1 — n^/n cq ) leading to Eq. 0 is only valid if ( is in kinetic equilibrium, or at least 







if /c(pO//c(p) = ~ p ')/ T . It is possible for the goldstino to elastically scatter off of the 

thermal bath as given by the interactions in Eq. (|3]): if the elastic interaction rate is large 
enough, the goldstinos would be in kinetic equilibrium with the thermal bath. However, 
because the £ + ■£’( A —>■ £ + i/s(A M ) process is suppressed by F k while the singlc-goldstino 
production channels are suppressed only by iy, kinetic decoupling takes place before chem¬ 
ical decoupling]^] To be precise, the momentum distribution of £ is determined not only 
by elastic scattering but also by chemical interactions. Thus if the production happens at 
T rh, energy-momentum conservation just tells us that the momentum distribution of 
thermally produced £ would be the form of /£(p) oc e~ p ^ T . Also if the chemical interactions 
are efficient enough (thus, for small F£), the produced goldstinos will still have a equilibrium 
distribution function, /£(p) = /f q = exp(— p/T). These arguments are not sufficient to 
justify the form of equation around the time of decoupling, since as we will see, at low T 
with small iy-term, goldstinos with different momentum decouple at different temperatures 
and, in the absence of elastic scattering, do not re-thermalize. Furthermore, the continuous 
entropy injection during matter dominated era causes the goldstinos decoupled earlier to 
be colder than those decoupled later. Therefore it is necessary to solve the non-integrated 
version of the Boltzmann equation for the distribution function /^(p): 


Substituting Hdt = din a, for T < fh this can be rewritten as 


( 12 ) 


as c a/ c 


d In a d In p 


= 1 


k \ t \ 


i-p/T 


Hp 2 J 


exp 


P_ 

T 


ml 

1H-^ 

4 p 2 


+ (£>—)■ A, "0 —t Ap) 


(13) 


in the limit of m^ —> 0. In the Appendix, we provide the Boltzmann equation for non- 
negligible m$ < rh. In the following, we will only consider decays from one generation of 
squarks, q, with mass rriq of O(TeV) while the others are assumed to be heavier. Adding 
other contributions is straightforward. For example, production by gaugino decays will take 
the same form with the substitution rriq —> nig and changing the number of sfermions to 
number of gauginos. For simplicity, the number of massless degrees of freedom (?*(T) is 
taken constant, g* = 85, in the whole range T R <T < rriq. 

3 This is the opposite behavior than for WIMPs, in which even after chemical decoupling WIMPs elastically 
scatter off of the thermal bath and remain in kinetic equilibrium until lower temperatures. 
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The freeze-out temperature for T/. a < rriq is calculated from 



M P TI 


f.o. 


1 gen. 


< r ^C?7 %nf/ 



( 14 ) 


where = m|/(167rF^), and the result is 



(15) 



This value is no longer true for large F-term. In this case, the freeze-out temperature 
becomes well above rriq, and it is mostly determined by scattering process. The freeze-out 
abundance is quite diluted by entropy production, and freeze-in production dominates the 
goldstino abundance as shown in Fig. [lj 

In the following, we derive the goldstino yield for each case of freeze-in and out. 

B. Freeze-in 

First, we will consider the simpler case of intermediate or high F^, for which goldstinos 
never reach chemical equilibrium for T < rriq ; their abundance is gradually increased by the 
thermal decay process (ii) until it is not efficient, after which they are diluted during the 
rest of the matter-dominated era. 

For n ^ <C n®- q , the resulting yield at reheating can be computed as: 



15.6 Mp iy> c . 


L 


(16) 


where sr = (27r 2 g*/45)Tft, and tji is the time at a = cir. K a is the modified Bessel function 
of the second kind, and Ki(x)/K 2 {x) ~ 1 — 3/(2a;) for x 1. tj is the initial time and 


we took it at the freeze-out of chemical interactions. The integrand on the second line 


shows high powers of (T R /T) caused by entropy injection and temperature dependence of 
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the freeze-in production rate, and because the Boltzmann suppression factor at low T from 


n~ q , the production rate is most efficient around T ~ rriq/ 10. The third line of Eq. (16) is 


obtained assuming T R -C m g -/10 -C Tf „_. Given the expression (15) for Tf ai this corresponds 
to Fq (500 TeV) 2 for rriq = 1 TeV. For lower values of F^, the freeze-out abundance is 
more important, which will be treated next. Here we just note that this yield is sizable, and 


will discuss late-time implications in Section IV 


C. Freeze-out 


For smaller /y, one could solve the same Boltzmann equation, (J9|)-(p^Tj) , and find 


tiq \ ec i 


^ ) = 1.7 x 10" 6 I — 




FO 


85 


-R 


lOGeV 


5 / 50GeV x 5 


Tf. 0 . 


(17) 


This would be incorrect, because the RHS of the Boltzmann equation ( [TTj ) was found assum¬ 
ing that the goldstinos are in kinetic equilibrium with the rest of the thermal bath. When 
goldstinos are injected in the bath via decays of non-relativistic particles, their momentum 
distribution is peaked around rriq/2. Because the 2 4 2 elastic scattering is frozen out at 
a higher temperature, it does not thermalize the distribution function. The result is that 
goldstinos with high momentum easily inverse decay back into superpartners, while goldsti¬ 
nos at low momentum are effectively frozen out. Therefore one can expect the correct value 
would be smaller than Eq. ©• 

This behavior can be understood explicitly by looking at the RHS of the Boltzmann 


equation for the distribution function, Eq. (13). We can estimate the effective ratio between 


the production rate and the expansion rate as 


ip, a) = 


Tq^ q mqT 2 


1 gen. 


Hp 3 


exp 


ATp 


(18) 


For a given scale factor a (corresponding to a given temperature T(a)), R^(p,a) changes 
with the momentum; in particular, if R^{p,ao) <C 1, the goldstinos with momentum p are 
decoupled, while if R^(p, a 0 ) 3> 1 they are in equilibrium. In Fig. [2j we plot R^ as a function 
of p for different temperatures, fixing rriq = 1 TeV, = (100 TeV) 2 , T R = 20 GeV. For 
example, at T = rriq/5 = 200 GeV all goldstinos with momentum p > 50 GeV are in thermal 
equilibrium, while at T = mg/10 = 100 GeV only the goldstinos with momentum p > 
150 GeV = 1.5T are, with the goldstino whose momentum is smaller than the temperature 
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FIG. 2. Left: Effective ratio between the production rate and the Hubble parameter, as a function 
of the goldstino momentum and for different temperatures. Interactions are in chemical equilibrium 
when 1, while the low-momentum region is frozen-out. Right: Ratio between the goldstino 

distribution function and the equilibrium distribution function, at rriq/T = 30. The continuous 
line comes from the numerical integration of the Boltzmann equation, while the magenta dashed 


line is the approximate analytical solution of Eq. (22). 


all decoupled. Thus at T « rriq/10 the result is an earlier departure of the number density 


from its equilibrium value compared to that of assuming kinetic equilibrium, Eq. (15) (this 
earlier departure can also be seen in Fig. [6] in the Appendix). 

Because Rq becomes very small once it is below one, we can solve the Boltzmann equation 
with zero RHS at low momentum, 


dfc(p,a) _ df ( (p,a ) 
d In a d In p 


for p < Pf. 0 .(a), (19) 


with boundary condition f^(p,a) = e p /A a ) f or p > p^ o ( a ). Here Pf . 0 .( a ) i s th e freeze-out 
momentum at a given temperature T(a), defined by R^{pf_ 0 _(a),a) = 1: 


P/.o.O) 


(k m g ) 2 

T(a) 


( 20 ) 


1/2 


= 0.13 


62.3 


62.3+ 4 In 


50 GeV 


+ 8 In 


(100 TeV) 2 


+ 2 111 ^ + 

9 * 


; In 


Tr 


20 GeV 


+ 24 In AA 
K c 


( 21 ) 
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Finally, the solution to the Boltzmann equation is 


fdP,a) = 


1 exp 

- (p/r c (“)) 8/11 

( exp[ 

-P/T(a% 


V > Vf.o{a ) 


T<(a)=(— Y /3 (—), 

W mqj \ a ) 


( 22 ) 

(23) 


To illustrate this better, in Fig. [2] we show the ratio between the resulting goldstino 
distribution function and the equilibrium distribution function, oc e _P//T , for a given tem¬ 
perature T = rriq /30 ~ 33 GeV. In blue, we show the results from numerical integration of 
the Boltzmann equation; in dashed, the analytical expression for the distribution function, 


Eq. (22), is shown, in good agreement with the numerical results. It is seen that the gold- 
stinos at high momentum are in equilibrium, while the low-momentum ones are suppressed. 
The late-time yield at low temperatures T < Tr, is given by 


^ — 6.8 x 10 -7 (— 


S / FO 


\9* 


’R 


10 GeV 


130 GeV 

hi^rriq 


(24) 


This is shown in Fig. [3j where we also show the result following from the kinetic equilibrium 


assumption, Eq. (17), and the abundance found without considering the matter-dominated 
epoch, that is the case Tr = Tmaxi hi which the Boltzmann suppression of the superpartner 
number density results in a negligible goldstino abundance at low reheating temperatures. 
Comparing the blue and magenta lines, we can conclude that a naive treatment of the 
Boltzmann equation overestimates the abundance by a factor of about 3. 

As a reminder, these results were found in the limit of small rri^, but they do not change 
much if the goldstino mass is sizeable. Even for rri( = 100 GeV, the final yield changes only 
by about 10%. In the Appendix, we show the full Boltzmann equation for the massive case, 
as well as numerical results for m^ up to 200 GeV (see Fig. [6]). 

Finally, as the distribution function deviates from the kinetic equilibrium case, one can 
also consider if the goldstinos produced would form colder or warmer dark matter, when 
compared to the case in which the particles are in kinetic equilibrium. The average momen¬ 
tum for the goldstino at T < Tr can be evaluated as 


(p) ~ 26 T c (a) = 0.36 T 


9*(T) \ 
9* ) 


1/3 


T, 


R 


10 GeV 


5/3 


130 GeV 


R^Ulq ) 


5/3 


(25) 


This should be compared to the thermal averaged value (p)t — 3 T. For low reheating 
temperatures, goldstinos are colder than the background temperature. This is due to two 
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FIG. 3. Late-time goldstino yield as a function of the reheating temperature T R , given = 
(100 TeV) 2 and rriq = 1 TeV. The blue line represents the numerical results, which coincide with 


the analytical result of Eq. (24), while the red line is the result that one would have found if kinetic 


equilibrium was assumed. The dashed red line is the yield found neglecting the matter-dominated 


era. 


competing effects: at first, they are produced at a higher momentum, p ~ mg/2, after which 
they are redshifted between production and reheating. For low T R , the second effect is 
dominant. 


D. Non-Thermal Production 

Goldstinos can also be produced non-thermally, for example by direct moduli/inflaton 
decays, or by squark (or other lightest WIMP particles in the MSSM) decays after freeze- 
out, for which n^/s = riq/s. The former is a model-dependent effect and it can be sizable 
or not. For what concerns the latter, it can be important for large F-terms: first of all the 
life-time of squark (with rriq = 1 TeV) has to be short enough to decay before BBN, implying 
F( (10 5 TeV) 2 [£]. Then, if > (5 x 10 4 TeV) 2 , the life-time of squark is long enough to 
decay after squark freeze-out by pair annihilation. In this case, the resulting energy density 
of the goldstinos is given by the non-thermal contribution 

rFh 2 = — -Qgh 2 . (26) 

rriq 

Since the reheating temperature is lower than the freeze-out temperature of squark anni¬ 
hilation, TjP*, Qgh 2 is also diluted by a factor of (Tr/TJ/*) 3 compared to usual freeze-out 
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abundance with Tr —> oo [2]. Because of the substantial model-dependence in the non- 


thermal goldstino abundance, the results in Eqs. (16)-(24) should be considered as conserva¬ 


tive results, as it is always possible to produce more goldstinos by introducing non-thermal 
processes. 

IV. LATE-TIME IMPLICATIONS 


The goldstinos produced in the early matter-dominated era are generally lighter than 
any other superpartner, except the gravitino. As such, they can provide a meta-stable dark 
matter candidate, even for very low reheating temperatures, Tr ~ 1 GeV. The lifetime for 
the decay to a gravitino, ( —> was computed in Ref. [B] as: 

F c \ 2 /100 GeV N 


7£ ~ 10 22 sec 


(27) 


(100 TeV) 2 / \ rriQ 

Although this is typically larger than the age of the universe, to — 10 17 sec, indirect detection 
limits on decaying dark matter are more stringent, with lower limits > C>(10 26 —10 27 ) sec 
for dark matter decaying to quark-antiquark pairs [THH2U] . As a crude estimate, we take the 
same order of magnitude, r™ n ~ 10 26 sec for the limits on the decaying goldstino. These 
can be avoided with goldstini lighter than 100 GeV, or larger F<y 


A. Assuming F-parity 

If F-parity is conserved, the goldstino is effectively stable in most of the parameter space. 
For small F-terrn (corresponding to goldstinos produced before freeze-out), the present dark 
matter density is 




IFO 


= 0.19 


\ f 85 


- R 


130 GeV 


V. 


(28) 


. 1 MeV J \g*J VlOGeVy \ k^rriq J 
The allowed region is shown in Fig. [4j and spans the range 0.5 GeV < Tr < 30 GeV. The 
result is only logarithmically dependent on the increase of F^, until Fq = (500 TeV) 2 , for 
which Tf 0 > niq/ 10. For larger F-terrn, the relevant process is freeze-in, and the dark 


matter abundance can be evaluated from Eq. (|16|): 

(Sj c /! 2 ) F /= 0.11 ( m< 


( 8 1\ 


3/2 


2MeV7 \gj V 10GeV 


Ti 


/ (500 TeV) 5 

V F C 


1 TeV 


nin 


(29) 
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FIG. 4. Late-time energy density of goldstinos, with small F-term Fq = (100 TeV) 2 and mq = 
1 TeV, for different values of the goldstino mass. For mq < 1 keV the goldstino is warm or hot 


dark matter. For mq in the range 100 GeV — 1 TeV (yellow region) the goldstino lifetime, Eq. (27), 
is too short, r < 10 26 sec, and is excluded by DM indirect detection constraints. The horizontal 
black line marks the observed value of the DM abundance, Rdm^ 2 = 0.12. 


To summarize this section, the observed dark matter abundance can easily be produced 
in a matter-dominated era, at temperatures above a low reheating temperature Tr. 


B. Assuming F-parity violation 

Another well-motivated possibility is that of F-parity violation. In the general case, 
baryon and lepton number would be violated and the proton would be unstable. Neverthe¬ 
less, if lepton or baryon number were to be independently conserved on their own, proton 
stability would be achieved accidentally. In the following, we will discuss the case in which 
baryon number is violated while lepton number is conserved (the other case will have a 
similar phenomenology). This is also interesting because baryonic F-parity violation can 
account for the matter-antimatter asymmetry of the universe with low reheating tempera¬ 
tures, while other scenarios (such as leptogenesis) require higher temperatures. The baryonic 
RPV operator in the superpotential is 

Wbrpv = -^u\d)d c k + h.c ., ( 30 ) 
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where the contraction of the color indices with an e abc tensor is understood and as a conse¬ 
quence j 7 ^ k. 

In low-scale gauge mediation with a single SUSY breaking sector, the gravitino is very 
light and there are still proton decay channels of the type p —* iT0 3 / 2 , mediated by RPV 
interactions. The resulting limits are very stringent and were discussed in [211 [22 ]. In the 
case of multiple goldstini, the proton could potentially decay to any goldstino lighter than 
1 GeV. This is particularly dangerous when is small, independently of the goldstino 
mass: for example, the weakest limit is 


A^g < 1.31 x 10 


-6 


m r , 


-X 


c 


(31) 


.ITeV/ (lOOTeV ) 2 

To avoid this constraint, we require that rriQ > m p \ however, now it is the goldstino which 
is unstable, as the decay channel ( —>■ Uidjdk is open. The lifetime is 

, in3 / i i / 10 GeV x 9 ' - w x 2 

tq =1.57 x lCr sec 


A ijk , 


rriQ 


( nig \U F ( 

V ITeV / \ (100 TeV ) 2 


(32) 


where A"- fc is the largest RPV coupling for which the decay is kinematically accessible. As 
the goldstino mass naturally lies in the interval 1 — 100 GeV, the top quark is not accessible 
and the most relevant operator with few constraints from flavor physics is A^UYs' 


The lifetime (32) of the goldstino naturally falls in a range that is probed by Big Bang 
Nucleosynthesis: if a large amount of energy is injected during the thermal plasma during 
BBN, the primordial abundance of light elements is changed and would go against observa¬ 
tions. In particular, the case of hadronic decays was studied in great details in Refs. [24] [27] . 
In the following, we will use the results of Ref. [55], where limits on the abundance M x Yx 
of a decaying particle X were set in the lifetime range 1CT 2 sec < tx < 10 12 sec, for different 
masses M x = 100 GeV, 1 TeV, 10 TeV. 

As we are also interested in particles with lighter masses, we need to extrapolate their 
results to M x = 10 GeV and below. Therefore, we will shortly review the source of the 
limits. For short lifetimes (r < 10 2 sec), the mesons and nucleons produced by X thermalize 
quickly, and the main consequence of the decay is the increase of the neutron-to-proton ratio, 
n/p , resulting in larger abundances of D and 4 He. For longer lifetimes, mesons decay and 
primary protons and neutrons scatter inelastically off of the background nuclei, generating 
hadronic showers, dissociating 4 He and producing D, T, 3 He, which also result in higher 
amounts of 6 Li, 'Li. At r > 10 7 sec, the neutrons decay away and only protons are left, with 
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a smaller effect on 4 He-dissociation. On the other hand, electromagnetic decay products 
( 7 , e + , e ~) are thermalized by processes such as 7 + 7 bg —» e + + e _ if their energy is above 
the threshold E t h = m 2 /22T [25] : one should compare this threshold to the binding energy 
of D and 4 He, respectively 2.2 and 28.3 MeV: if it is higher, non-thermalized photons will 
dissociate deuterium and helium. As the bath temperature decreases with time, photo¬ 
dissociation of D and 4 He is active for r > 10 4 sec and r > 10 6 sec, respectively. 


We simulate the goldstino decay with Pythia 8.2 [29] and get the total number of charged 
particles and EM energy per Q decay, for different values of m<y We then translate the results 
of Ref. [25] to lower masses. For the sake of simplicity, we only use the dominant constraints, 
that is primordial helium abundance ( Y p ), deuterium to hydrogen ratio (D/H) and helium-4 
to deuterium ratio ( 4 He/D). Our results are shown in Fig. [HJ where for comparison we also 
show the constraints of Ref. B3- 


For most of the lifetime range, the limits on are of order 10 12 — 10 14 GeV. Com¬ 


paring this to the yields found in Eqs. (16) and (24), only very low reheating temperatures 
are allowed. In Fig. |5j we show contours of the maximum allowed reheating temperatures 
in the rn^ — A"- fc plane: apart from a small shaded region in the top right corner, for which 
77 < 10 ” 2 sec, the upper limit on the reheating temperature is of order 0.5 GeV. In the lower 
left corner, we shaded the area in which 77 ^ 10 12 sec and goldstino decays do not interfere 
with BBN: here they would be constrained by diffuse X-ray and 7 -ray emissions, and the 
limits on m^Y^ are typically more strict than the BBN ones. As this discussion is beyond 
the scope of the present paper, we only display the BBN limits as a conservative bound. 


In the presence of A-parity violation, there is one more channel for the production of gold- 
stini that had not been analyzed so far in the literature: the scattering qq —> q( does not 
require any on-shell superpartners, in contrast with the A-parity conserving processes dis¬ 
cussed above. Using the interaction in Eq. Q, the freeze-out temperature for the scattering 

is 


rri RPV 
1 f.o. 


70 GeV 


V85/ 


(lOOTeV ) 2 


2/5 / mg \4/s / 1 \ 2 / 5 /10 GeV\ 2/5 

VTbV/ vxj V t r J 


(33) 


The goldstino distribution is the same as that in kinetic equilibrium as long as they are in 
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FIG. 5. Left: Upper limits on the energy density of a hadronically decaying particle X , with, 
from bottom to top, Mx = 10 GeV, 100 GeV, 1 TeV, 10 TeV, evaluated from Ref. [251 ■ Right: 
Maximum reheating temperature allowed by BBN constraints while varying X'-- k and ?n<p The 
other parameters are fixed as = (100 TeV) 2 and rriq = 1 TeV. 


chemical equilibrium, and once we solve the relevant Boltzmann equation we find the yield, 


7 J<\ 
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= 7.0 x 10 -7 |A"| 2 In 


RPV 


rnRPV 

1 fo- 

Tr 


85 A 

<]*) 


3/2 


T, 


r y/ (100TeV) 2 y /TeVA 


/ 85 A 3/2 / T r A 7 /(100 TeV) 2 V / 
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\ m q J 

4 


V m q 


(34) 


Here the hrst line shows the production during the matter-dominated era and the second 
line is the production after reheating; in general, the former dominates over the latter, and 


the T-parity conserving freeze-out contribution from Eq. (24) is larger than both. For the 
RPV production rate to be sizable, the only option would be Tr = Tmax •C m, for which 
there is no early matter-dominated era. In this case, only the last line in the equation 
above contributes to goldstino production. There are still strong limits from BBN, and 
the maximum reheating temperature is of order 1 — 10 GeV, with smaller RPV couplings 
allowing slightly larger reheating temperatures {T^ ax ~ 20 GeV for X" = 10~ 5 ). 

The limits on the reheating temperature cited so far corresponded to F ,[ = (100 TeV) 2 , 
with goldstinos generated at freeze-out. The dependence of T^ ax on larger T-terms is hrst 
logarithmic, until F\ [ ~ (500 TeV) 2 , and then scales as T^ ax oc (T^/(500 TeV) 2 ) 2 / 7 when the 
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freeze-in contribution (16) becomes dominant. For example, = (10 6 TeV) 2 corresponds 


to Tft ax ~ 100 GeV for = 10 GeV and \” jk = 1. 

We end this section by summarizing the strong bounds on the reheating temperature 
in the case in which R-parity is violated. For the freeze-out case, the maximum reheating 
temperature allowed by BBN constraints is of order 1 GeV. 

Some implications for what concerns baryogenesis can be drawn: because any baryon 
asymmetry produced at higher temperatures in the matter-dominated era will be diluted 
away, the baryon asymmetry should be generated at temperatures between T r and Trbn ~ 
10 MeV. This is possible in the LSP baryogenesis scenario of Ref. [30] if the goldstino decays 
before BBN (in the upper right corner of Fig. [5]): for example, with parameters chosen as 


F< = (100 TeV) 2 , A» t ~ 0.1, the goldstino abundance in Eq. ([24]) can be large enough for 
baryogenesis if ~ T R ~ 50 GeV. 

Even though we have taken rriq = 1 TeV as a benchmark point, it is worth noting that 
current LHC limits on RPV squarks are less stringent: for a light top squark decaying 
to three quarks, the CMS collaboration excludes masses up to 350-385 GeV in Ref. [31] • 
Gluino limits are stronger: in Ref. [32], CMS excludes gluinos in the decay channel g —> tbs 
up to 900 GeV, while in the same channel ATLAS excludes gluinos up to 874 GeV, with 
limits of order 800 GeV for different flavor composition of the final states. Thus, squarks 
are still allowed to be lighter than in our benchmark point. 


V. CONCLUSIONS 

In this work, we have discussed thermal production of goldstinos as an example of super- 
weakly interacting particles during an early matter-dominated era ending at reheating. This 
is important when the reheating temperature is low, in the GeV range or below, as particle 
production through an heavier sector (superparticles) occurs only at higher temperatures 
that were not achieved during radiation domination. 

We have analyzed in detail the production of an uneaten goldstino by solving the Boltz¬ 
mann equation for its momentum distribution function, and revisited the cosmological impli¬ 
cations. When the goldstino is stable enough, thermal production can provides the correct 
dark matter density even for reheating temperatures as small as 1 GeV. If R-parity is vi¬ 
olated, the goldstino has to be heavier than the proton and is meta-stable, with a lifetime 
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range that naturally interferes with BBN. In this case, reheating temperatures higher than 
1 GeV are excluded for almost all of the small parameter space. 

Such low reheating temperatures suggest a low scale of inflation and/or introducing cer¬ 
tain symmetries to prevent coupling between the inflaton and visible fields. For example, if 
the inflaton decay rate to the MSSM is Planck-suppressed, Tj = m 3 /Mj,, the reheating tem¬ 
perature is Tp ~ 1 GeV(m//10 3 TeV) 3//2 . In this case we can also find an upper bound on 
Tmax from its definition in Eq. (J7]), T MAX < 10 3 TeV(Tp/ GeV) 2//3 . Since the non-thermal 
production is proportional to Tfi/rrij oc T// 3 , it could be more important at low reheating 
temperature, and full analysis considering a specific inflation model is needed. In this pa¬ 
per, we presented thermal production of goldstino as model independent contributions. Our 
results are conservative bounds on the abundance, because the non-thermal productions of 
SWIMP are just additive quantities. 

We showed that the goldstino is one of many good examples in which a momentum 
dependent process gives a sizable difference compared to that assuming thermal distribution 
of the momentum, even if they are produced from thermal bath. Such effects are also 
discussed in Ref. [S3] for the production of sterile neutrino dark matter, and in Ref. |33j for 
leptogenesis from heavy Majorana neutrino decays. Since the period of kinetic decoupling 
and production mechanism are also important to the evolution of density perturbation of 
dark matter, the study of the full Boltzmann equations in the perturbed spacetime can give 
observable consequences for small scale structures. 
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Appendix A: Boltzmann Equation for Massive Goldstino 

In this section we derive the Boltzmann equation for the momentum distribution func¬ 
tion of (, _//(p), keeping the dependence on a non-zero mass m q. The dominant source of 
production is sfermion decays, (j) —>■</ + Vh when the temperature is lower that the sfermion 
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mass fhy,, and the contribution from elastic scattering is ignored. As we did in the rest of this 
work, we will consider the Boltzmann equation and distribution functions in the classical 
limit, i.e. assuming / x ( p) < 1, and / x q (p) — e~ Ex ^ T , where E x = sjrn^ + p 2 for particle x- 
The corresponding Boltzmann equation in the limit of massless quarks is 
dfc_df c N , n d ft_ nn 

Hp W ” l/d 


9<t>9 y 
~2E~c 


(Sj%(^)%; (2 ’ r)li<1,(p * -p* -p") - uu), (ai) 

where p = |p|. Tree level T symmetry ensures \A4 t p^.^\ 2 = \A4 ( ^ ( j > \ 2 at leading order. 
Sfermions are in thermal equilibrium, which is maintained by the interactions with the 
background SM fields (0 + 0* yy + A(f> + 'ijj —> 0 + 0, ...) before such interactions 
are frozen. After that, the distribution of sfermions is determined by the interaction with 
the goldstinos, and we need to solve the coupled Boltzmann equations. However, most of 
goldstinos are produced before the freeze-out of sfermion-SM interactions, so we can safely 
take /^(Py) = e~ E ^l T . 

Using the identity /^ q 5 ( 4 ) (j^ - - p l ) = /^ q /^ (4) (pj - to represent after 

integrating over p,/, , we get 


k 



dpcj) „eq P (j> 


E, 

P\ 


W* 


I d cos 9^6 (E# 

g-E^/T 


E, 


C~Pp) 


(A2) 


I dpfy/d cos 6,% | ' 

where p^ = sjp^ + p 2 - 2 p^p cos 9^, so |<40/dcos600 = p^p/p^- For given p, the integral 
domain D\p\ is determined by E^ — E^ = \Jp 2 ( f > + P 2 — 2p<f>P cos 9^ for —1 < cos 9^ 

Then, we hnd p^ < p<p < p^, where 

oA E cTP) ( E c ±p) 


< 1 . 


Py = m: 


2 m 2 


(A3) 


In terms of energy variable E^, E^ < E^ < E+, where 


IS = + AlM, 


In massless limit, —> 0, we obtain 


2 m 2 
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FIG. 6 . Evolution of the goldstino number density for non-negligible m^. Here we have fixed 
m.q = 1 TeV, = (100 TeV ) 2 and Tr =40 GeV. From top to bottom, the continuous lines are 
the yield computed assuming kinetic equilibrium in the massless case, and the yields using the 
full Boltzmann equation, for three different values of the masses, m £ = 0 GeV, 100 GeV, 200 GeV. 
The dashed lines are equilibrium number densities for the same masses. 


which was used in Eq. (13). It is straightforward to evaluate C[fc\ as 

c[/d = 


" /T 


n 1 

Finally, the Boltzmann equation can be written as 

a/c % 


k \ 


\e~ E ~J T e ~K/ T ] 

e~ E d T ) 

V 9(E ( p ) 
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d In a d In p 
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( V w q mqT\ 

—E7/T -E+/T 
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This expression replaces Eq. (13) in the case of non-negligible m £. 

For example, in Figure [ 6 ] we show the evolution of the goldstino number density for 
different values of the goldstino mass for a fixed T rt = 40 GeV. From top to bottom, with 
the continuous lines we show the yield computed assuming kinetic equilibrium and = 0 , 
and the yields computed solving the above Boltzmann equation for three different values of 
the masses, = 0 GeV, 100 GeV, 200 GeV. We find that the final yield of goldstinos is not 
changed much for rriQ < 100 GeV. The dashed lines are the equilibrium number densities for 
the three different masses: we see that the effect of the masses is much smaller than what 
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naively expected by looking at the equilibrium number density. 
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